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ABSTRACT 

We have detected low-amphtude radial-velocity variations in two stars, USNO-Bl.O 
1219-0005465 (GSC 02265-00107 = WASP-1) and USNO-Bl.O 0964-0543604 (GSC 
00522-01199 = WASP-2). Both stars were identified as being likely host stars of tran- 
siting exoplanets in the 2004 SuperWASP wide-field transit survey. Using the newly- 
commissioned radial-velocity spectrograph SOPHIE at the Observatoire de Haute- 
Provence, we found that both objects exhibit reflex orbital radial-velocity variations 
with amplitudes characteristic of planetary-mass companions and in-phase with the 
photometric orbits. Line-bisector studies rule out faint blended binaries as the cause 
of either the radial-velocity variations or the transits. We perform preliminary spectral 
analyses of the host stars, which together with their radial-velocity variations and fits 
to the transit light curves, yield estimates of the planetary masses and radii. WASP-lb 
and WASP-2b have orbital periods of 2.52 and 2.15 days respectively. Given mass esti- 
mates for their F7V and KIV primaries we derive planet masses 0.80 to 0.98 and 0.81 
to 0.95 times that of Jupiter respectively. WASP-lb appears to have an inflated radius 
of at least 1.33 i?jup, whereas WASP-2b has a radius in the range 0.65 to 1.26 i?jup- 

Key words: methods: data analysis - stars: planetary systems - techniques: radial 
velocities - techniques: photometric 
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1 INTRODUCTION 



Extra-solar planets that transit their parent stars are of 
key interest because their masses and radii can be deter- 
mined directly, providing; clues to their internal compo- 
sitions l|Guillot et al.ll2006h . They define the mass-radius- 
separ ation relation for irradiated giant planets (|Mazeh et al.l 
I2OO5I ). They provide unique insig hts into their thermal prop- 
Charbonneau et all [2OO5I : iDeming et all l2005l : 



erties (e.g. 



iDeming et al.l 2006 ) and the chernical compositions of their 
atmospheres ([Charbonneau et al.ll2002l : IVidal-Madiar et al.l 
l2003l : IVidal-Madiar et al.ll2004l ). 



Table 1. Target list for SOPHIE radial- velocity foUowup pro- 
gramme. The ISWASP identifiers give the J2000 stellar coor- 
dinates. The second column gives the number N^i of SOPHIE 
observations that were required to classify and eliminate each 
non-planetary system, and the total number of observations for 
planet-bearing systems. Half the astrophysical false positives were 
eliminated after a single observation. Five more were eliminated 
after a second observation, leaving 8 narrow, single-lined targets. 



ISWASP ID 



SOPHIE foUowup classification 



The first exoplanet found to exhibit transits , HD 
209458b (|Charbonneau et al.|[200ol : iHenrv et aLll2000l ). was 
initially discovered using the radial- velocity method that has 
to date yielded the vast m ajority of the 210 known exo- 
planets. Its inflated radius (iBrown et all I2OOII ) presents a 
challenge to theories of the structure and evolution of irra- 
diated exoplanets. New radial-velocity discoveries are rou- 
tinely subjected to careful photometric foUowup at around 
the predicted times of transit. Indeed, recent radial- velocity 
surveys targeting bright stars of high metallicity with the 
specific goal of discovering new hot Jupiters ha ve revealed 
two new transiting plane ts in the last two years (jSato et al-l 
l2005l : lBouchv et ahlboosh . Transiting planets thus comprise 
roughly 10% of all the hot Jupiters with orbital periods un- 
der 5 days or so, which is consistent with expectations for 
randomly-oriented orbits. 

The complementary approach is to look for transits 
first, then to seek evidence of their planetary nature in 
radial-velocity foUowup observations. The OGLE project pi- 
oneered this approach, which has to date yielded five transit- 
ing exoplanet candidates for which radial- velocity variations 
have been detected with amplitud es ind icative of planetary- 
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mass companions (iKonacki et ah ' 2003; B ouchv et al. II2OO4 



iKonacki et"alll2004l : IPont et al.ll2004 : .Konacki et al.ll2005^ 



Several teams have embarked on ultra-wide field 
searches for transiting exoplanets. In this complementary 
approach, small-aperture CCD imaging systems coupled to 
commercial camera optics secure light curves of millions 
of stars. The TrES survey has s o far yielded two transit- 
ing planets (lAlonso et all |2004| : lO'Donovan et al.1 l2006"bl ). 
The XQ survey and HAT surveys have so far found one 
each (jMcCullough et al.l I2OO6I : iBakos et al.l[2006h . Here we 
present the first two transiting exoplanets from the Super- 
WASP survey. The hardware, data analysis pipeline and 
ar chive methodology for this project are described in detail 
bv lPoUacco et al.l (|2006l ). 



In a series o f rec ent papers, IChristian et~aLl l|2006l ). 
IStreet et aP ()2006t ) and lLister et al.l ( 20061 ') presented transit 
candidates from the inaugural 2004 May-September observ- 
ing season of the five SuperWASP wide-field survey cameras 
that were operating at that time. The light curves of some 
1.1 X 10^ stars in the magnitude range 8 < V < 13 were 
searched for periodic shallow transits, and several dozen 
plausible candidates were identified for detailed foUowup. 
The candidate selectio n methodology is descr i bed in detail 
in these papers and bv lCoUier Cameron et al] (|2006l ). 



2 SOPHIE AND OFXB OBSERVATIONS 

We conducted a radial- velocity survey of a sample of high- 
priority SuperWASP transit candidat es, using the newly - 
commissioned SOPHIE spectrograph (jBouchv et al.l I2OO6I ) 
on the 1.93-m telescope at the Observatoire de Haute- 
Provence during the four nights from 2006 August 31 to 
2006 September 3, and between September 12 and 19 dur- 
ing the science verification phase of SOPHIE. SOPHIE is a 
bench-mounted, fibre-fed spectrograph bui lt on the same de - 
sign principles as the HARPS instrument (|Pepe et al.ll200j ) 
on the ESO 3.6-m telescope at La Silla. The spectrograph's 
thermal environment is carefully controlled, with the aim of 
achieving radial- velocity measurements with stability better 
than 2 m s^^. For the targets studied here, a radial-velocity 
precision of 10 to 15 m s~^ is adequate to establish or re- 
ject the planetary nature of a transit candidate. We there- 
fore elected to use SOPHIE's High-Efficiency (HE) mode, 
which has resolving power A/ A A = 35000. The CCD de- 
tector records 39 spectral orders spanning the wavelength 
range from 387 to 694 nm. Radial velocities are determined 
by cross-correlation with a mask spectruin matched to the 
spectral type of the target l|Baranne et al.|[l996l : IPepe et al.l 
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Table 2. Journal of radial-velocity measurements of WASP-1 and WASP-2. The ISWASP identifiers give the J2000 stellar coordinates 
of the photometric apertures; the USNO-Bl.O number denotes the star for which the radial-velocity measurements were secured. The 
uncertainties given here include 10.0 m s~^ systematic error in added in quadrature to the formal photon-noise error. The fifth and 
sixth columns give the FWHM of the CCF dip and the contrast of the dip as a fraction of the weighted mean continuum level. The 
signal-to- noise ratio near 550 nm is given in column 6, and the spectral type of the cross-correlation mask in column 7. 

HJD icxp Vr FWHM Contrast S:N Mask Notes 

(s) km s"-' km s"-*- % Sp. type 



ISWASP J002040. 07-1-315923. 7: USNO-Bl.O 1219-0005465 = GSC 02265-00107 = WASP-1 
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I2OO2I ). Automatic data reduction at the telescope allows 
highly efficient candidate selection and data assessment. 

The target l ist wa s drawn from the ca ndid ate papers of 
IChristian etaP (|2006l l. lStreet et all (|2006l l and lLister et all 
l|2006ll , and from further candidate lists covering different 
regions of the sky for which candidate papers are currently 
in preparation. Table [T] gives a brief summary of all targets 
observed and the outcomes of the SOPHIE observations. 

A full investigation of the natures of the other objects 
observed during this radial-velocity study is currently in 
progress. The results of this investigation and its implica- 
tions for improving our pre-selection criteria will be pub- 
lished in a companion paper. Here we present the two targets 
in the survey sample that have single, narrow-lined cross- 
correlation functions (CCF) with radial velocity variations 
in agreement with an oscillation precisely phased with the 
ephemeris predicted from SuperWASP data. The amplitudes 
of these radial velocity oscillations are less than a few hun- 
dred m s~^ and no significant line-bisector variations are 
detected (see Section [5]). Thus we can conclude that each 
of these two stars harbours a transiting planet. The radial- 
velocity measurements for the two stars are listed in Table 
1. 

A complete transit of WASP-2 was observed using a 
CCD camera with iZ-band filters on the 60-cm telescope 
of the Observatoire Frangois-Xavier Bagnoud at St-Luc 
(OFXB), on the night of 2006 September 12/13 UT. A full 
transit of WASP-1 was observed on the morning of 2006 
October 2 UT using an SBIG STIOXME CCD camera with 
7?-band filter on the 0.35-m Schmidt-Cassegrain Telescope 
at the Volunteer Observatory at Knoxville, Tennessee. 



3 STELLAR PARAMETERS 

The extracted SOPHIE spectra were used for a preliminary 
analysis using the uclsyn spectral synthesis package and 
ATLAS9 models without convective overshooting (Castelli 
et al. 1997). The Ha, Nal D and Mgl b lines were used 
as diagnostics of both Tbs and logg. The abundances do 
not appear to be substantially different from solar. We used 
these values to infer the radii and masses of the stars, as 
listed in Table |31 Fo r WASP-1, compariso n with the stellar 
evolution models of ICirardi et al. I l|2000t ) gives maximum- 
likelihood values M* = I.ISM© and R* = 1.24iZQ, but 
many models with 1.06 < M,/Mq < 1.39 and 1.04 < 
Rt/RQ < 1.92 satisfy the spectroscopic constraints, because 
the main sequence is very wide in that temperature range. 
The radius and mass estimates for WASP-2 are better con- 
strained. We have also used the available BV and 2MASS 
photometry to estiniate Te S using the Infrared Flux Method 
(|Blackwell fc Shallis|[l977l ). which gave results in agreement 
with that obtained from the spectral analysis. 



4 PLANETARY PARAMETERS 

The formal precision of the radial- velocity observations de- 
pends on the signal-to-noise ratio of the spectrum, the sharp- 
ness and density of the stellar lines, and the scattered- light 
background in the instrument. The formal errors on the ve- 
locity measures are given by the semi-empirical estimator 
o-Rv = 1.7\/FWHM/(S:N*Contrast). The Contrast param- 
eter quantifies the contrast of the CCF peak against the 
additional background signal from light leakage into the as- 
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Table 3. Stellar parameters derived from preliminary spectral 
analyses with U CLSYN. Masse s and radii are derived from the 
Padua models of lGirardi et alj 1I2OOOI) assuming [Fe/H]= 0.1 ±0.2 
to reflect the uncertainty in the metallicity. 



Parameter 



WASP-1 



WASP-2 



GSC 

WASP V (mag) 
Spectral type 
Tcff (K) 
logs 

My (mag) 

M,/Mq 



02265-00107 00522-01199 



11.79 
F7V 
6200 ± 200 
4.3 ±0.3 
3.9 ± 0.4 

1 24+"-^* 



11.98 

KIV 
5200 ± 200 
4.3 ±0.3 
6.2 ± 0.5 
n 70+0.15 

0.78 ±0.06 



yet incomplete spectrograph enclosure. Because we did not 
use simultaneous thorium-argon wavelength calibration, the 
accuracy of the radial-velocity measurements is limited by 
the stability of the spectrograph over the 2 to 3 hours be- 
tween successive thorium-argon calibration exposures, con- 
sidering that its thermal control was not yet optimized. Tests 
performed during the commissioning of the instrument indi- 
cate that the velocity drift during a night is typically in the 
range ±10 m s~^. Taking into consideration additional un- 
certainties coming from the wavelength solution of the H.E. 
mode and the guiding noise, we estimated that during our 
run, the systematic RV errors were 10 m/s. Although the 
wavelength calibration may vary only slowly during a given 
night, the drifts may reasonably be expected to be uncor- 
related from one night to the next. Guiding noise should 
be uncorrelated even between successive observations. Since 
most targets were observed only once per night, we treat 
the additional systematic error as uncorrelated. We there- 
fore used the quadrature sum of the formal and systematic 
errors (Table El for all model fits. 

The transits of WASP-1 and WASP-2 can be timed with 
a precision of about 20 minutes from the 2004 SuperWASP 
data set. The first and last transits in this data set are sep- 
arated by about 120 days, so the accumulated uncertainties 
in the transit timings at the epoch of the 2006 SOPHIE, 
OFXB and Volunteer Observatory observations are no more 
than a few hours. There is no ambiguity in number of cycles 
between the 2004 and 2006 data sets. The OFXB and Volun- 
teer Observatory transit observations establish the improved 
photometric ephemerides in Table |4] 

We estimated the stellar and planetary radii and the 
planetary masses by minimising for the photometric and 
radial-velocity measurem ents simultaneously wi th respect 
to the analytic model of iMandel fc Agol ()2002l ) for small 
planets {Rp/R, < 0.1), assuming that the planets have cir- 
cular orbits. In modelling the SuperWASP photometry we 
used linear limb d arkening coefficients it — 0.51 and 0.63 
Ivan Hammelll993l ) for WASP-1 and WASP-2 respectively. 
For the more precise data from OFXB and Volun t eer O b- 
servatory we used the 4-coefficient model of IClarej (|2000"). 

We used a bootstrap analysis to estimate the errors on 
the parameters of the model. This entails taking the residu- 



Table 4. Results of simultaneous minimum-x^ circular-orbit fits 
to the photometric and radial- velocity data for WASP-1 and 
WASP-2. The parameters of the lightcurve model are given in 
terms of the radius of the star and planet (ij* and Rp, respec- 
tively) , the separation of the stars (a) and the inclination (i). The 
total number of degrees of freedom, including the Af^y radial- 
velocity measurements, is N^{. The contribution of the radial- 
velocity data to the value of is Xrv' ^^ii^S the errors given 
in Table [2] Data in transit are given reduced weight. Standard 
errors on the parameters are derived using a bootstrap analysis 
as described in the text. 





WASP-lb 


WASP-2b 


Transit HJD 


2453912.514 ± 0.001 


2453991.5146 ±0.0044 


Period (days) 


2.51995 ± 0.00001 
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als from the optimum light-curve fit, applying an arbitrary 
phase shift and restoring the model transit at phase zero. 
This preserves both outliers and the correlated noise char- 
acteristics of the WASP data. The synthetic radial-velocity 
data are generated by sampling the best fit radial-velocity 
curve at the observed phases and adding gaussian random 
deviates with the appropriate standard error. These syn- 
thetic data are then fitted repeatedly, to recover the distri- 
butions of the fitted parameter values. 

There is a strong degeneracy between the impact pa- 
rameter b — acosi/i?* and the parameter R,/a when fitting 
planet transit lightcurves of the quality presented here. We 
therefore present the results of minimum-^'^ fits with the 
value of b fixed at values b = (corresponding to the mini- 
mum radius of the star) and b — 0.8 for both stars. Models 
with b > 0.8 give noticeably worse fits to the lightcurves. 

The radial- velocity measurements and minimum-^'^ fits 
are shown as a function of orbital phase in Fig. [T] The sinu- 
soidal variation in radial velocity is clearly seen. The Ught 
curves and minimum-x'^ fits to the phases around the tran- 
sits are shown in Fig. E) The parameters of the best fits 
are given in Table |4l An additional 50 m s~^ was added in 
quadrature to the uncertainties of radial-velocity observa- 
tions during transits , to allow for the R ossiter-McLaughlin 
effect jRossite dligli : lMcLaughlinlll924j ) . 

In spite of this, the contribution of the radial-velocity 
data to the total value of in both stars is larger than 
expected. WASP-1 yields a spectroscopic Xa = 11-6 for 5 
degrees of freedom, while WASP-2 gives Xs ~ 13-4 for 7 
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Figure 1. Radial velocities of WASP-land WASP-2 plotted against barycentric Julian date (upper) and folded on the photometric 
ephemerides (lower). In all panels the solid line represents the best-fitting circular-orbit solution. 



degrees of freedom. More extensive observations will be re- 
quired to establish the cause of this small additional "jit- 
ter" in the radial-velocity measurements. A common cause 
of such radial- velocity "jitter" is rotationally-modulated dis- 
tortion of the line profiles arising from magnetic activity in 
the stellar chromosphere and photosphere. While the S:N in 
the bluest orders of the SOPHIE spectra is insufficient for us 
to measure the chromospheric emission cores in the Ca II HK 
lines in any meaningful way, both stars have narrow CCFs 
that suggest wsini values less than ~ 5 km s~^ and hence 
low to moderate activity levels. We allow for this additional 
variability in the bootstrap error analysis by adding 12 m 
s~^ of additional radial-velocity jitter to the synthetic radial- 
velocity data. We note that this amount of jitter is within 
the ranges determined empirically for F and K stars with low 
to moderate levels of chromospheric activity f Wrightl l2005h . 
and conclude that chromospheric activity is a likely cause. 

Also given in Table|4]are the masses and radii of WASP- 
Ib and WASP-2b derived using the best-fit model parame- 
ters. Despite the ambiguities in the light-curve solution it is 
clear that the data presented show that WASP-1 and WASP- 
2 have planetary-mass companions with gas-giant radii. 



5 FAINT BLENDED-BINARY SCENARIOS 

Although many types of stellar binary system can mimic a 
planet-like transit signal, most are easily eliminated. Strong 
rotational broadening implies tidal synchronisation by a 
massive companion. Grazing double or single-lined stel- 
lar binaries reveal their nature after one or two observa- 
tions. Triple systems can, however, mimi c transiting planets 
in a way that is difficu lt to eliminate l|Torres et al.l |2004| : 
lO'Donovan et al.l [2006al ). A bright single F, G or K star 
with a much fainter, physically associated eclipsing-binary 
companion will apparently produce shallow eclipses. If the 



WASP-1 



1.04 - 




0.96 - 



-0.10 -0.05 0.00 0.05 0.10 

Photometric Phase 
WASP-2 



1 .04 - 




0.96 - 



-0.10 -0.05 0.00 0.05 0.10 

Photometric Phase 

Figure 2. Transit profiles of WASP-1 (upper) and WASP-2 
(lower), fitted with the optimum limb-darkened models with pa- 
rameters as in Table |4] Small symbols denote 2004 season Super- 
WASP photometry; filled circles show the Volunteer Observatory 
light curve of 2006 October 1 (WASP-1) and the OFXB light 
curve of 2006 September 12 (WASP-2). 



bright primary has narrow absorption lines, the faint bi- 
nary's lines will produce apparent periodic velocity shifts in 
the combined spectrum, in phase with the photometric or- 
bit. The apparent shift occurs because the broadened lines 
of the tidally-synchronised primary of the eclipsing binary 
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are Doppler shifted by orbital motions into the wings of the 
composite line profile. 

Several tests can reveal this type of system. Because the 
spectral types of the bright star and the cooler eclipsing bi- 
nary are very different, cross-correlation masks of different 
spectr al types tend to yie ld different orbital velocity ampli- 
tudes ijSantos et al.ll2002l l. We computed the CCF with dif- 
ferent masks without significant change in the radial- velocity 
values. 

Because the apparent Doppler shift arises through a 
variable asymmetry in the line wings, line-bisector analysis 
of the cross-correlati on function is also a n effective probe for 
this type of system (jOueloz et al.ll200ll ). When the fainter, 
moving spectrum is broadened by rotation, the induced ve- 
locity shift is sufficiently weak that the apparent velocity 
variation tends to be small unless the binary's lines are 
strong enough to give a clearly-visible asymmetry in the 
combined line profile. The change in line-bisector velocity 
from the wings to the core of the line tends to be greater than 
the induced velocity amplitude. We measured the asymme- 
tries of the cross-co rrelation func t ion p eaks using the line- 
bisector method of lOueloz ^~el\ (j200ll ). For WASP-1 and 
WASP-2 we found the scatter in bisector velocities from the 
wings to the core of the CCF profile to be substantially 
less than the measured orbital velocity amplitude, and un- 
correlated with orbital phase. A third candidate, ISWASP 
J2103f8.01+080117.8, failed this test and was eliminated as 
a planet candidate. 

Therefore, we can eliminate the scenario "bright single 
star plus faint, late-type short-period eclipsing binary" with 
confidence for WASP-1 and WASP-2. While other more in- 
tricate scenarios are possible, we were not able to contrive 
any that could explain both the photometric and the veloc- 
ity signals while remaining credible. 

At the 380 and 140 pc distances of WASP-1 and WASP- 
2, it should be possible to resolve such binary companions 
at separations of a few tens of AU or more with the help 
of adaptive optics. As an additional line of defence against 
this type of astrophysical false positive, we secured high- 
resolution H-hanA images of both targets with the NAOMI 
adaptive-optics system on the 4.2-m William Herschel Tele- 
scope on the nights of 2006 September 6 and 7. Images 
with corrected FWHM = 0.25 arcsec reveal that WASP-1 
has a stellar companion 4.7 arcsec to the north and 3.7 
magnitudes fainter at H . The SOPHIE fibre aperture has a 
diameter of 3 arcsec. Seeing and guiding errors may allow 
some of the companion's light into the fibre, but we can be 
fairly confident that the contamination is not significant at 
visible wavelengths. Even if the companion were an eclips- 
ing binary, it would be unlikely to mimic the radial- velocity 
signature of a planet, and indeed the line- bisector analysis 
eliminates this possibility. NAOMI H-band images taken on 
2006 September 7 with corrected FWHM=0.2 arcsec, and 
using the OSCA coronagraph system, show that WASP-2 
has a stellar companion 2.7 magnitudes fainter at H located 
0.7" to the east. This falls within the SOPHIE fibre aper- 
ture. Additional NAOMI images secured during transit on 
2006 September 10 20:00 to 20:20 UT with 0.2 arcsec cor- 
rected FWHM showed no sign of the ~1.5-mag deep eclipse 
in the companion that would be needed to mimic a transit. 
Future AO observations should reveal whether these faint 
stellar companions (which are the only objects visible be- 
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Figure 3. Mass-orbital separation diagram for the 14 known 
transiting exoplanets. WASP-lb and WASP-2b follow the general 
trend with only high-mass planets surviving at small separations. 
(From http://obswww.unige.ch/^pont/TRANSITS.htm) 



sides WASP-1 and WASP-2 in their respective NAOMI fields 
of view) are chance alignments or common-proper-motion 
companions. 



6 DISCUSSION AND CONCLUSIONS 

We have detected the presence of radial-velocity varia- 
tions in two exoplanetary transit candidates. Our prelim- 
inary analysis yields masses between 0.8 and 1.0 Mjup 
for both planets. WASP-lb and WASP-2b lie between 
the "hot Jupiters" and the "very hot Jupiters" in the 
mass-orbital separation diagram for transiting hot Jupiters, 
(Fig. [3]). Their intermediate masses appear consistent with 
the general trend to ward high mas s es at the smallest orbital 
separations noted bv lMazeh et all ()2005l ). WASP-2b in par- 
ticular lies close to the minimum separation at which planets 
in this mass range survive, making it a good candidate for 
future mass-loss studies. The radius of WASP-lb, which or- 
bits an F7V star, is poorly constrained by the SuperWASP 
data alone, but appears to be at least 1.33 i?jup. WASP-lb 
seems thus to be an expanded, l ow-density planet, similar 
to HD 209458b and HAT-P-lb teak os et al] l2006l l. Addi- 
tional photometry of WASP-2b, which orbits a Kl dwarf 
at a slightly smaller orbital separation yields a radius close 
to that of Jupiter, suggesting a substantially higher density. 
Additional high- precision photometry is needed to refine the 
radii and densities of both planets; indeed, precise radius es- 
timates hav e recently been derived from high-precision pho- 
tomet ry by IShporer et al.l l|2006l) and ICharbonneau et ahl 
(l2006h . confirming the oversized nature of WASP-lb. 
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